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When addressing the thermodynamics of finite-sized systems, one must specify whether one 
wants to fix conserved charges to a sharp value or whether one is content to fix their thermodynamic 
average. In other words, contrary to the thermodynamic limit, different statistical ensembles are not 
equivalent. When treating the plasma phases of gauge field theories perturbatively in the canonical 
ensemble, unexpected new difficulties arise in comparison with the usual grand canonical treatment. 
The purpose of this paper is to expose these difficulties and show how they can be remedied, thus 
(^-^ , recovering a well-defined description of plasmas even in a finite-size, canonical ensemble setting. 

■ For definiteness, a specific model is considered, namely QCD1+1 with SU(2) color; however, the 

^\ ' treatment presented should be applicable also to higher-dimensional systems with different types of 

Coulomb interaction. 
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PACS number(s): ll.10.Wx, 11.15.-q, 12.38.Cy, 12.38.Mh 



I. INTRODUCTION 



CN . 

The fundamental interactions governing microscopic behavior in nature are described by gauge theories. These 
i theories generically give rise to long-range Coulomb interactions, and for this reason plasma phases represent an 
almost universal phenomenon in dense systems, whether one is considering quarks and gluons in a heavy ion collision 
or the electron gas in a metal. Many of the plasmas accessible to experiment are quite small, not only in the case 
of the heavy-ion collision |l]], but also e.g. in mesoscopic systems ||], which have become a wealthy source of new 
phenomena in the past years. In such a setting, finite-size effects become important, and they are indeed one of the 
causes of the high diversity of mesoscopic physics. 

The most obvious effect of a finite configuration space is that the corresponding momentum space becomes discrete. 
However, there is also a slightly more subtle source of new behavior connected with the ensemble one uses to describe 
the statistical properties of the system at hand. Commonly, in the presence of a conserved charge N, one uses the 
^H. grand canonical ensemble to calculate the partition function, 

^ I Z GC = Tr e-KB-i**) (1) 

>■ 

T *~J 1 where fi is the corresponding chemical potential, (3 the inverse temperature, and H the Hamiltonian. One may at the 
end of the calculation fix /i to impose a desired mean value of N . Physically, this corresponds to the system being 
able to exchange energy and particles with its surroundings. In many cases, however, this is an idealization, and the 
system is only able to adjust its charge to a limited extent, or not at all. The prime example is the colliding heavy-ion 
system, which must remain an exact color singlet at all times. In such a situation, the canonical ensemble is more 
appropriate, 

Z C = [ daTr e -W-*"Vf-"o)/f>) (2 ) 



This restricts the trace to states in which the charge N takes exactly the value No- In the following, No will always 
be taken to be zero, without loss of generality. 

Whereas for very large systems, (j|) and @) give identical results because the fluctuations of N in the grand canonical 
ensemble die out as the inverse root of the size, in a finite system ([!]) and (|J) are different. Now, at first glance one 
may think that the finite size corrections are easy to obtain if one has, as usual, calculated (|l]); simply put in an 
imaginary chemical potential (i = ia/(3 and do an integration to obtain (Q). Such considerations have been made 
for ideal gases e.g. in ||- ||. A closer look however reveals that, at imaginary chemical potential, a plasma does 
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not conform to the recipes conventionally used in its formal treatment. Specifically, the venerated prescription of [Q, 
namely summing the so-called ring diagrams to a logarithm, fails because one runs into the cut of the logarithm. The 
main formal advance of the present treatment is to show how this difficulty can be remedied. Thus it will be possible 
to regain a well-defined description of the plasma even in a finite-size, canonical ensemble setting. 

The specific model to be considered here is QCD in one space dimension in the high temperature or high baryon 
density perturbative regime S. It must be emphasized however that the formal developments presented are quite 
general and are straightforwardly applied also to higher-dimensional systems with any type of Coulombic interaction. 
A one-dimensional model is chosen because in this setting, the infrared divergences appear most clearly, uncluttered 
by the usual ultraviolet problems. The results however are expected to enjoy a wide range of applicability. The goal of 
the treatment is to evaluate the contribution of the plasmon effect to the thermodynamics of a finite system, exhibiting 
especially the differences in behavior away from the thermodynamic limit depending on which ensemble is used. Of 
course, the size of the system is still taken to be large enough to permit a physical interpretation as a plasma, i.e. if g 
is the coupling constant in the relevant Lagrangian and L is the length of the system, then gL should still be large^J. 
Otherwise, the concept of screening the Coulomb potential becomes meaningless and the perturbative expansion of 
the theory is in powers of g 2 L 2 . Physically, large gL means that the plasmons, whose size is controlled by 1/g in 
one space dimension, still fit well into the length L. Correspondingly, in three space dimensions, the plasmon size 
behaves (9) like 1/gT and thus the relevant dimensionless parameter which must remain sizeable would be gTL. Only 
under such conditions is it meaningful to consider the usual plasmon contribution to the thermodynamic potential 
and evaluate its behavior as 1/gL is increased away from zero. 

II. THE MODEL 

The axial gauge Hamiltonian of QCDi + i with two colors reads in momentum space 

9 2 f dq. a/ 



H = J dp X Hp)(p l5 + m 10 )x(p) + ^.f(q)3 a (-q) (3) 

where the fermion fields x carry Dirac and color indices and the SU(2) currents j a are 

j a {q)= J d PX Hp)\x{p + q) (4) 
with the Pauli matrices a a . For later use it is also useful to introduce the U(l) currents 

f(q) = \J d PX \p)x{p + q) (5) 

The Hamiltonian (^|) conserves the U(l) charge j°(0), corresponding to the baryon number, and the SU(2) color charges 
j a (0). For definiteness, it will be assumed that the plasma may exchange baryon number with its surroundings, but 
not color charge. This will illustrate all the possible complications; all alternative scenarios can be treated along 
similar lines. In particular, it is straightforward e.g. to specialize to QED in the canonical ensemble. Thus, the goal 
is to evaluate the partition function 

Z= I dG{a) Tr exp [~/3(H - M j°(0) - ia a f(0)/p)] (6) 

JSU(2) 

where dG{a) denotes the Haar measure of SU(2) [^)- [||. Generically, a further simplification is possible at this point. 
The trace is invariant under a global unitary transformation of the fermion fields, 

X -> X' = U X (7) 

and also the Hamiltonian and the baryon number remain unchanged under this transformation (this is a global residual 
gauge invariance still present after gauge fixing). Therefore, one may choose to diagonalize a a cr a , 
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U^a a a a U = ffVK) 2 + (a 2 ) 2 + (« 3 ) 2 = ^ (8) 
i.e. the partition function depends only on the length of a, not its direction in color space, 

z = h r da sin2 ? Tr exp ^~ p{H ~ w '° (o) ~ wj3(o)//?) ] = i / * ^ sin2 1 z(a) (9) 

In general, only the (imaginary) chemical potentials associated with the Cartan subalgebra of the gauge group are 
relevant. In addition to the symmetry exploited above, there still remains a discrete translational symmetry in a and 
also a parity symmetry corresponding to the freedom of permuting the eigenvalues in the diagonalized matrix in (^). 

In the following, ([)]) will be referred to as the "canonical" partition function (even though the U(l) charge is treated 
grand canonically) . On the other hand, treating also the SU(2) color charges grand canonically (this will be referred 
to as the "grand canonical" case) corresponds simply to dropping the average over a in (g) and setting Z = Z(a = 0). 

It should be mentioned that in the specific setting of QCD1+1, the Hamiltonian (||) is only valid in the color singlet 
sector. In general, the Hamiltonian of QCDi + i contains additional couplings to residual gauge degrees of freedom Jl0| , 
the dynamics of which eliminate all non-singlet states from the spectrum. Thus, the partition function of QCDi + i 
strictly speaking contains only singlet states already by virtue of the dynamics. Of course it is fully equivalent to 
use the color singlet Hamiltonian (Jfy and impose the color singlet constraint by hand as in (||). On the other hand, 
for purposes of comparison below also the grand canonical partition function corresponding to the Hamiltonian (|^) 
will be considered. Strictly speaking, this is not a physically meaningful alternative. Therefore, from the point 
of view of QCDi + i what is being compared is a physically correct calculation (namely, the canonical one) with the 
idealized grand canonical one which is often adopted in the literature due to technical convenience. Of course, in other 
settings, such as electrodynamical plasmas in three space dimensions, both alternatives may be physically meaningful, 
depending on the laboratory conditions. 

III. IDEAL GAS 

Momentarily neglecting the Coulomb interaction, one can now immediately read off the ideal gas partition function 
\nZ (a) = A J dpJ2 [Ml + e-^+^+^/W) + ln(l + e -/3(«»-M/a-*of«a/V)jj (10) 

i 

= ^Jdp [ln(l + 2cos(a/2) e ^ + ^ 2 )+e- 2/3 ^/ 2 )) 

+ ln(l + 2 cos(a/2)e-' 3(£ ^ /2) + e -20(«*-l«/2))l (11) 

where e p = \Jp 2 + m 2 . In some special cases, this can be evaluated explicitely as long as one still assumes the 
momenta to be continuous. E.g. in the classical nonrelativistic limit m 3> m — fi/2 ^> T one obtains 



lnZ (a) = ^cos(a/2) J^e-^^/ 2 > (12) 

7T V P 



whereas in the limit m — fi — one has p2[ 



7r a 2 \ 



ln ^)=^--— ) (13) 

Taking now into account the interaction perturbatively, assume that lnZ(a) has been calculated to some order in 
g (note that one usually calculates directly hiZj(a) by considering only connected Feynman diagrams), 

\a.Z{a) =]nZ (a) +lnZ I (a) (14) 

According to (||), this must now be averaged over a, 

lnZ = bx— [ da sin 2 -Z (a)Zi(a) (15) 
2tt J 2 

, 1 f, ■ 2 a r 7 r x I da sin 2 (a/2)Z (a)Z I (a) 

= ln— / da sin -Z a (a)+ln — — (16) 

2ir J I J da sin (a/2)Z (a) 
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Thus Zo(a), apart from giving the ideal gas contribution in the first term of (|16[), also acts as a measure for averaging 
Zi(a). Now, inspection of ([l2| ) and (|l^), or in general ([Til), reveals that Zq(ol) falls off as a Gaussian with a width of 
l/\/L in a for large L. Therefore, in the thermodynamic limit, only the vicinity of a = contributes to InZ. Thus 
one recovers the well-known result that the canonical and grand canonical ensembles are equivalent for large systems, 
or, formulated in another way, color singlet constraints become irrelevant in the thermodynamic limit. 

For finite L, on the other hand, the two ensembles differ. Whereas the grand canonical ideal gas partition function 
is obtained simply by setting a = in (|Tlj), ([l2|), or (|T3|), in the canonical ensemble one must average over a. In 
special cases, this can be accomplished analytically. In the classical nonrelativistic limit, the average over ( |l2| ) can be 
carried out to give 



\nZ = \n(2I l (x)/x) with x = -J^L e -^ m -^ (17) 

7T V P 



where I\{x) is a Bessel function of imaginary argument fl^| . In the case m = /i = one obtains using (13) 



lnZ„ = g+ln 



' ' rf(2vrVi) - e~& (Re crf(27rV^ + r-^) - Rc erf(-^) 



(18) 



where erf denotes the error function and x — L/(AnP). The comparison between the ensembles is illustrated for the 
ideal gas in figures (jl]) and (j^). Also, to give an idea how the effect of discretizing the momenta compares with the 
effect of varying the ensemble, the free partition function has been calculated using antiperiodic boundary conditions 
for the fermion fields, i.e. 

(2n+ 1)tt f , 2tt 

as well as periodic boundary conditions, 



P^Pn = —, d P^Y^ ( } 

Evidently both figures, though stemming from vastly different regions of the phase diagram, display a dominance of the 
effect due to varying choice of ensemble except at very low L. Formally, what happens is that the difference between an 
integral and the corresponding Riemann sum behaves as l/£, as can be inferred from the Euler-McLaurin summation 
formula |r|. Thus, merely discretizing the momenta leads to an expansion of the form InZo ~ O(L) + 0(1) + . . . for 
the logarithm of the partition function. On the other hand, the averaging procedure involved in doing a canonical 
calculation in general introduces power corrections in L into the partition function Zq, as can be explicitely observed 



in ( |l7| ) and (18). This leads to an expansion of the form InZo ~ O(L) + O(lnL) + . . . upon taking the logarithm. 
At low L, the question which effect dominates depends on the details of the dynamics, as evidenced in figures (0) 
and (||). Whereas in figure (j|), the discretization effects are negligible even at LT = 1, in the ultrarelativistic case 
the discretization effects become comparable to the effect of varying the ensemble at sufficiently low L. The effect of 
using a canonical ensemble is reinforced by using antiperiodic boundary conditions and weakened by using periodic 
ones. A similar comparison will be possible for the plasmon contribution calculated in the next section. 



IV. PLASMON CONTRIBUTION 



This section is concerned with the main object of the present treatment, the perturbative evaluation of the plasmon 
contribution to the thermodynamic potential. Whereas this effect arises only at order 0(g 3 ) in gauge theories in three 
space dimensions jxifl , it already contributes at order O(g) in one space dimension due to the different volume element 
in momentum space. Thus it represents the dominant perturbative effect here. 

As a side remark, note that in QCDi + i, this is in fact the only term of the perturbative expansion which can be 
obtained in terms of Feynman diagrams. Already in the next order of perturbation theory, the Feynman expansion 
breaks down due to irremediable ambiguities introduced by the infrared divergences |Q. Formally, one can give 
convergent sums over selected subclasses of diagrams such that virtually arbitrary fractional powers (greater than 
one) of the coupling constant g arc generated. On the other hand, it is known e.g. in the classical nonrelativistic 
limit of QCDi + i |L5| that expanding the exact equation of state in the coupling constant yields a series made up 
exclusively of integer powers of g. Similar infrared problems arise in QCD 3+ i at order 0(g 6 ) jlq] . In the case of QED, 
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the infrared divergences are not as severe and the Feynman expansion can be carried out at least to some higher order 
in g. 

The perturbative expansion of the expression (|^) for the partition function may be obtained by standard methods 
jnj [|l7|. The corresponding Feynman rules can be read off from (JoJ) and (|J) as follows: 



• Fermion propagator with momentum p, color i and Matsubara frequency oj s — (2s + 1)tt//3: 

i uj s - ifi/2 + af i a/2[3 + i-f 5 p — i-f m 



Ki(p, s) 



(3 (lu s - in/2 + <4a/2/3) 2 + p 2 + m 2 



(21) 



Coulomb interaction with momentum transfer q and energy transfer 2nr//3 (cf. figure (|J)): 



y p \ -* 

16tt 



(22) 



All internal momenta, frequencies, color and Dirac indices must be integrated or summed over, respectively; fermion 
loops are associated with an additional factor (—1), and diagrams with N Coulomb interactions have an overall factor 
1/NL Every connected part of a diagram in addition receives a factor L/2ir; note that these rules correspond to 
already having used momentum conservation at all vertices. By considering only all connected diagrams, one directly 
obtains the logarithm of the partition function, an extensive quantity in the infinite volume limit. 

The plasmon contribution IiiZr is obtained by summing the ring diagrams (cf. figure (||)). Physically, this is 
interpreted as a screening of the Coulomb interaction by the production of fermion- antifcrmion pairs; thus the infrared 
divergences resulting from the long-range character of the interaction are dampened. Since there are 2 N ~ 1 (N — 1)! 
ways of combining N vertices to a ring, the Feynman rules give 



lnZ R (a) = — £ 



dq Tr 



color 1 ( n 0. 9 = 0) 



^ N 

N=l 



N 



(23) 



where r is the integer labeling the frequency flow around the ring, q is the momentum flow around the ring, and the 
polarization II is a matrix in color space, 



n aa (r,0) 



47T/3 



1,3 



E / d p 



[(<. 



2f3 u ti> 



2/3 a jj > + £ : 



(24) 



Note that in (|2^), only the value of the polarization at q — is taken. Higher orders in the Taylor expansion of II 
around q — dampen the infrared divergence and thus ultimately lead only to terms of higher order in the perturbative 
expansion in g 111 ] . The sum over s in (E4) is readily carried out by standard methods [|17| , 



n aa (r, 0) 



Gji{a) 



9 ^°**«*r -f(^-4) 



32tt 



31 



dp 



1 



cosh § (e p + f + igfoj.) cosh|(e p + f + i&of.) 



'2P"33> 2 ' 

1 



cosh §(e p 



^ - i— a 3 ) cosh S.(e — £ 

2 L 213 a ]3> Whn 2 2 



(25) 



7 — a 3 ) 

L 20° ill 



(26) 



Note G22 = Gti an< i ^21 — G12 > 0. Furthermore, it is advantageous to diagonalize II in color space, since 
for the trace over arbitrary powers of II. One arrives at the eigenvalues 



calls 



n x (r) = 2 5 2 ReGn(a) 



sin ■ 



sin — 

n 2 (r) = 2 5 2 G 12 (a) 2 - 



n 3 (r) = 2 5 2 G 12 (a) 2 - 



(27) 
(28) 
(29) 
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Why III is left in this form will become clear below. Note that in the case of electrodynamics, the only term which 
appears is of the form of IIi; the SU(2) case thus exhibits additional structures compared with the U(l) case. 

Cast as above, it becomes very transparent where the situation at finite imaginary chemical potential differs from 
the usual one. According to (E3h, one must now calculate 



J i r N=l vy ' 



N 



(30) 



Usually i.e. for a = 0, one would observe that (E7p-(E(J) only give contributions for r — and are negative there. 
Thus, the sum over r contains only one term, one may carry out the sum over N to obtain a logarithm, and finally 
the (now infrared- regular) q- integral may be done JtJ p7| . 

In the more general case considered here, it is not meaningful to sum to a logarithm, since IT2 and IT3 may have 
positive values depending on the value of r; also IIi, though only giving a contribution for r = 0, may become positive 
because ReGn turns to negative values at some finite a. However, as soon as the IIi become positive, one only has 
to go to sufficiently low q to run into the cut of the logarithm to which one may have naively summed ( |3(i| ) . 

There does not seem to be a physical interpretation for this failure. One is simply not doing a good job of defining 
what the initially ill-defined sum over N in (^) should be; the usual summation to a logarithm replaces one ill-defined 
quantity with another one. There is no physical principle which may e.g. tell one how to circumvent the cut in the 
logarithm by giving a small imaginary part to the momentum q. It must be emphasized that this problem has nothing 
to do with the question of how the infrared singularity in the Coulomb interaction (cf. (^)) should be defined. It 
already occurs for finite, sufficiently large L, where the momenta are discrete and where the Coulomb interaction 
consequently is completely well-defined (at least in the color-singlet sector considered here). The difficulty thus seems 
to be more related to the asymptotic nature of the perturbative expansion itself. 

However, one may indeed do a better job of defining (^) by proceeding in a different way, namely by exploiting 
the structure of (|2^)-(|29|) in r; this is the reason why Hi was not explicitely simplified to Hi = — 2S r o ReGn. The 
key observation is that 



E 



i 



(Trr - f ) N 



2 ^ |fc| f 



1 



(2ik) 



N 



(31) 



(this identity is derived in the appendix). In this way, one generates a crucial additional factor l/(N — 1)! which turns 
the sum over N in (30) into an exponential. Therefore one finally obtains 



hxZ R {a)- 



L x ^ 1 

Jki 

k^O 1 1 



dq e lke (e 



ike / Aikg sin(^)Re Gn(ct)/q 



1) 



+2e ika (e 4ik9 sin (f ) G i2( a )/« 



1) 



(32) 



This expression is completely well-defined. It represents the main improvement of the present treatment over more 
conventional ones. Note that the procedure above is reminiscent of the Borel summation method Jl8| , where one 
introduces a factor 1 = (A:!) -1 dt e~ l t k to generate better convergence of a series. The price one pays is the 
introduction of an additional integration. Here, on the other hand, the convergence-improving factor is generated 
simply by rearranging the already present r-summation. A certain price is implied by the poor numerical convergence 
of the fc-series as it stands. How this is best handled will be commented upon further below. 

In general, if the momenta are discrete, one must proceed from equation ( |32"| ) numerically. As long as the momenta 
are continuous, one may do the integral over q 12 (due to the reflection symmetry of ( p2[) in a and e around zero, 
these parameters will be taken positive in the following): 



In Z R (a) 



e tk \j\k\ sin - J Re Gn (a) I (-1 + isgn(fc) sgn(ReG n (a))) 



+2e""V|fc|sin-Gi 2 (a) (-1 + zsgn(fc)) 



(33) 



Furthermore, in the case a — * 0, one may use 

1 

m 



lim > — 1=& 



fc^o 



W^(-l + zsgn(fc)) = -2 r rf x Sin(X + 7r/4) =-2^F 
2 Jo V x 



(34) 



G 



and in this way verify that in a grand canonical treatment, or in a canonical treatment in the thermodynamic limit, 
one reproduces precisely the result of the usual prescription outlined directly after equation (|3^). As a side remark, 
note that taking additionally the classical nonrelativistic limit allows to evaluate G explicitely to give 

ln^ = -^L(^y /4 e -f(-^) (35) 

Combining this with the corresponding ideal gas result, determining p(/i,T) = dP/d/j.\T, and solving for /i to obtain 
the equation of state yields 

'-"•('-Tv^-) < 36 > 

This result was obtained in [l5| by a completely different method and thus provides a nice check on the present 
calculation. Similarly, one obtains in the case m = (still in the thermodynamic limit) 

^Z R = -\^ 9 L (37) 

Some further results for finite a are collected in the appendix. 

With the help of the general formula (|32|), one may now again perform a comparison of the different ensembles. 
Note that in practice, it is advantageous to perform the a-averaging over ([32] ) before doing the fc-sum. Then the 
fc-sum converges quite well above k ~ yL for large L, as can be inferred from stationary phase considerations. Of 
course, adopting this procedure forces one to treat rapidly oscillating a-integrands for large k, but this is numerically 
the lesser evil. Also, as before, allowance was made for a discretization of the momenta. Note that discretizing the 
momentum transfer q in ( |32[ ) , which can be interpreted as the plasmon momentum, always means taking 

27m f 2tt \ -> 

q -> q n = — , / dq -> — ( 38 ) 

regardless of whether one is considering antiperiodic or periodic boundary conditions for the fermion fields. Note 
also that the term q = is excluded from the sum. This is due to the fact that this term originates directly from 
the q = term in the original Coulomb interaction (cf. (g)), proportional to j a (0)j a {0)/q 2 - However, j a (0) is being 
constrained to be zero and therefore this term may be excluded from the Coulomb interaction. Physically, this simply 
takes account of the fact that for a neutral system, there is no physical content in a fluctuation translating all charges 
in the system uniformly. 

The comparison between the different cases is illustrated in figures (|^) and (^). The pictures are strikingly similar, 
even though taken from vastly different regions of the phase diagram. Again, at very high L, the effect of varying the 
ensemble dominates. However, in contradistinction to the ideal gas term, the effect of discretizing the momenta takes 
over already at quite high L. This does not come as a complete surprise. The simple argument used in the discussion 
of the ideal gas term, that the Euler-McLaurin series implies an expansion purely in powers of 1/L for the logarithm 
of the partition function if one discretizes the momenta, is not valid for the integral over the plasmon momentum q 
due to the singular nature of the integrand in ( p2] ) around q = 0. Thus, there is no simple argument anymore which 
would allow one to predict which effect dominates for large L. That the strong effect of momentum discretization 
here is indeed due to the plasmon momentum sum, as opposed to the fermion momenta contained in In Zq and the 
quantity G, is corroborated by the observation that figures (||) and (||) remain identical regardless of whether one 
uses periodic or antiperiodic boundary conditions for the fermions. Only at very low L, not displayed in the figures 
anymore, could one observe a dependence on the boundary conditions. This happens at such low L that the entire 
concept of treating the system as a plasma becomes invalid. Note that this is also true for the ideal gas contribution 
depicted in figures (|l|) and (|^). Below L/f3 w 10, it is impossible to simultaneously maintain gL 3> 1 and g/3 <C 1. 

Thus, the only relevant finite-size effects in the plasmon contribution stem from the use of a canonical ensemble and 
from discretizing the plasmon momentum. The two effects are roughly comparable in magnitude, with the ensemble 
effect dominating at very high L. 

V. DISCUSSION 

The present investigation has focussed on the different finite-size effects influencing the thermodynamics of a plasma. 
Special emphasis was placed on the evaluation of the plasmon contribution in a canonical ensemble as opposed to the 
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usual grand canonical treatment. Here, it turned out that the usual resummation prescription for the ring diagrams 
is too naive, and it was shown ho w the resulting ambiguity may be remedied via equation (|3l|). This latter trick, 



which in view of its derivation ( Al)-( A5) may be suggestively termed "Fourier regularization" , may be useful in quite 
diverse contexts as a means of properly defining ill-defined series. In the present context, it should be noted that 
the nontrivial manipulations purely involve the imaginary time direction or the corresponding frequency sums. Thus, 
while the present treatment focussed on a one-dimensional model, all manipulations are equally applicable for the 
analogous problem in an arbitrary number of space dimensions. Also, the versatility of the method is evidenced by 
its application to the more complicated SU(2) Coulomb interaction as opposed to the U(l) one. 

Using these formal developments, the partition function of QCDi + i for SU(2) color was calculated to order O(g), 
exhibiting the effect of varying the ensemble used and of discretizing the momenta. The dominant finite-size effects 
were identified to be the ones resulting from using a canonical ensemble and from discretizing the plasmon momentum. 
The effect of discretizing the fermion momenta in comparison was found to be negligible for interesting values of L. 
Of the two dominant effects mentioned above, the ensemble effect is strongest at very high L. 

Finally, it should be mentioned that there is an alternative interpretation of the imaginary chemical potential a over 
which one averages in the canonical ensemble: The configuration space of the model considered here is a cylinder, for 
finite L possibly closed to a torus depending on the spatial boundary conditions used. Moving on the cylinder parallel 
to its axis corresponds to moving in the space direction, moving around the cylinder corresponds to moving in the 
imaginary time direction. The parameter a can be interpreted as an (imaginary) background chromoelectric potential 
as would be produced by a hypothetical solenoid coinciding with the axis of the cylinder. This background potential 
introduces a phase into Green's functions which lead around the cylinder, or in more physical terms, the a-dependence 
corresponds to the Aharonov-Bohm effect induced by the solenoid on the thermodynamics. A similar situation was 
investigated recently in connection with the statistical properties of the Gross-Neveu model on a ring jl9| . One of 
the interesting results of this latter treatment was the observation of non-analytic behavior of the thermodynamical 
observables as a function of the Aharonov-Bohm phase. Disregarding for the moment that in the present calculation, 
the parameter a is in the end an integration variable, similar behavior is found here in quantities before a-averaging, 



e.g. in the first term in the square brackets of the expression (A7) for \n.Zu{a). It may be worthwhile to investigate 
whether there is physical meaning contained in this singularity. 
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APPENDIX: SOME USEFUL DERIVATIONS AND FORMULAS 



Identity (31) can be seen in the following way: 



I fdX"- 1 ^ 1 



E (7rr _|)iv {N -i)\\ da ) ^Trr-f (A1) 



N-l / A \ N - 1 



(N - 1)! \da 

2 n fj^N-lJ* 



(N-l)\ \da 



cot(a/2) (A2) 

oo 

^ sin ka (A3) 



fc=i 



2^ /d^ 1 



I^sgn(fc)e lfcQ (A4) 



(N-l)\ \daj 2i 

Note that the manipulations of Fourier series carried out above are well-defined in the context of the theory of 
distributions lEOfl. 
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As long as one considers continuous momenta, many quantities connected with the plasmon contribution to the 
partition function can be evaluated explicitely even for finite parameter a, in analogy with the ideal gas expressions 
given in equations (113) and (13). Thus, in the classical nonrelativistic limit one has 



G, 



SttV /? 



ReG n = cos(a/2)Gi 2 



(A6) 



and therefore, using (p3|), 



hxZ R {a) = — — 



1/4 



r (ro-/i/2) 



v /cos(a/2)6»(cos(a/2)) 



/° i~. — ", — ttttv^ sin(fca + 7r/4) 



(A7) 



where 9(x) denotes the step function. Note that the fc-sum can be expressed as a generalized zeta function [gT| . 
Sinilarly, in the limit m = fi = one obtains 



_ 1 a/2 1 

47rsm(a/2) An 



and from this 



In Z R (a) 



4a v sin(fca + 7r/4) 



7T ^ 



fc=l 



(A8) 



(A9) 
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FIG. 1. Ideal gas part of the logarithm of the partition function in the ultrarelativistic limit m = (X = 0. Short-dash-dotted: 
Grand canonical ensemble with discrete momenta (periodic boundary conditions); long-dash-dotted: Same, with antiperiodic 
boundary conditions. Solid line: Canonical ensemble with continuous momenta; short dashes and long dashes: Same with dis- 
cretized momenta, periodic and antiperiodic boundary conditions, respectively. A grand canonical calculation with continuous 
momenta gives a constant, namely the one all curves converge to at large L. 
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FIG. 2. Ideal gas part of the logarithm of the partition function in the classical nonrelativistic limit: m/T — 10, f/,/T = 14. 
Short-dash-dotted: Grand canonical ensemble with discrete momenta (periodic boundary conditions); long-dash-dotted: Same, 
with antiperiodic boundary conditions. Solid line: Canonical ensemble; in this case the curves obtained with continuous 
momenta and with momenta discretized in the various ways are indistinguishable. A grand canonical calculation with continuous 
momenta gives a constant, namely the one all curves converge to at large L. 



FIG. 3. Coulomb interaction 



FIG. 4. Sum of ring diagrams 

FIG. 5. Plasmon contribution to the logarithm of the partition function in the ultrarelativistic limit m = [i = 0, with 
g/T = 0.1. Dotted: Grand canonical ensemble with continuous momenta; solid line: Same, with discretized momenta. In 
the latter case, periodic and antiperiodic boundary conditions for the fermions give the same, shown, curve. Long dashes: 
Canonical ensemble with continuous momenta; short dashes: Same with discretized momenta. Again, the boundary conditions 
chosen for the fermions do not make a difference. 

FIG. 6. Plasmon contribution to the logarithm of the partition function in the classical nonrelativistic limit: 
m/T — 10,/i/T = 14, with g/T = 0.1. Dotted: Grand canonical ensemble with continuous momenta; solid line: Same, 
with discretized momenta. In the latter case, periodic and antiperiodic boundary conditions for the fermions give the same, 
shown, curve. Long dashes: Canonical ensemble with continuous momenta; short dashes: Same with discretized momenta. 
Again, the boundary conditions chosen for the fermions do not make a difference. 
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